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На основе применения преобразования простого растяжения-сжатия указана методика распространения решений задач, свя- 
занных с течениями несжимаемой вязкой жидкости в присутствии сферы, на варианты, когда сфера заменяется трехосным эл- 
липсоидом. Решена задача о медленном обтекании эллипсоида, указана простая расчетная формула для его сопротивления. По- 
казано удовлетворительное совпадение с литературными данными, соответствующими предельным случаям. 


Введение 

Некоторые задачи, связанные с интегрировани- 
ем уравнений Навье-Стокса 

Ѵ-у - ѵ х гоі ѵ = -Ей Ѵр - ^гоі гоі ѵ, сііѵ ѵ = 0, (1) 

или их упрощений, в случаях присутствия твердого 
или жидкого объекта в виде сферы единичного ра- 
диуса успешно решены путем использования функ- 
ции тока 

Ф=^ст 2 А(У)ып 2 Ѳ, ( 2 ) 

<Т =б,І+ЦІ + 1^к =<т[/С08Ѳ+8ІпѲ(УС08ф + Л8ІШД)]. 

При этом от функции Дет) требуется, чтобы она 
удовлетворяла уравнению 

гоі гоі гоі ѵ=гоІ гоі гоі гоі ( і//Ѵ ф) = 0, (3) 

что в общем случае влечет за собой представление 
/(<т) = с, +с 2 <т~’ + с 3 <у~ 3 +с 4 (т 2 . (4) 

При различном наборе констант интегрирова- 
ния в (4) функция тока (2) будет обеспечивать ки- 
нематические картины течения как во внутренних 
(вихри Адамара-Рябчинского-Хилла, с 2 =с 3 =0), так 
и во внешних (с 4 =0) областях. Для внутренних те- 
чений ускорение является консервативным векто- 
ром, и давление находится из полных уравнений 


движения Навье-Стокса; в идеальном внешнем по- 
токе (с 2 =0) давление определено интегралом Бер- 
нулли, а во внешнем вязком оно находится из урав- 
нений Стокса. 

Ниже обсуждается проблема распространения 
отмеченных и других решений на случай замены 
сферы трехосным эллипсоидом. Сопротивление 
эллипсоида вращения при его медленном движе- 
нии было найдено в [1] по теории ньютоновского 
потенциала притяжения. В [2] предпринята попыт- 
ка определить сопротивление набегающему потоку 
вязкой жидкости трехосного эллипсоида при па- 
раллельности потока и одной из полуосей путем 
использования преобразования простого растяже- 
ния-сжатия, при котором сфера переводится в эл- 
липсоид и наоборот. 

В общем случае трехосный эллипсоид имеет бес- 
конечно много сопротивлений, что зависит от его 
ориентации к набегающему потоку. Поэтому иссле- 
дование выполняется в предположении, что орты 
/ , у, к декартовой системы координат жестко связа- 
ны с главными центральными осями эллипсоида 



который в тексте будет фигурировать в виде урав- 
нения (при ( 7 = 1 ) 
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где I - линейный масштаб, используемый при 
обезразмеривании системы Навье-Стокса, а в ка- 
честве х, у, і участвуют теперь безразмерные коор- 
динаты. Используется линейный переход от раз- 
мерного радиуса-вектора г к безрамерному 

С- -■ г Г х - ѵ - 2 Г 

<т=$і +г\)+і;к=—і +— } + —к 

Л] А 2 


по схеме преобразования простого растяжения-сжатия 

г = хі + у] + гк = а—і +Ь~і+с —к = Ь + А 2 іу‘ + А^к ), 
а Ъ с 


= Лег, Л = 


(\ О 
О А, 

0 0 Аз] 


Особое внимание уделено нахождению реше- 
ния системы: 

Ѵр + гоі гоі ѵ = 0, ёіѵ ѵ = 0, гоі гоі гоі ѵ = 0 (5) 


в переменных вектора а в области 1<о<со при ус- 
ловии прилипания ѵ(1)=0. Необходимо отметить, 
что возможное решение в переменных г /, С, не яв- 
ляется универсальным, что можно пояснить следу- 
ющим примером. Линейное уравнение движения 
свидетельствует о том, что давление р является гар- 
монической функцией. Но, например, функция 


Р = ^ 


ѵ а 1 2 + а, 2 + а 2 




3 

2 


в пространстве і;, ц, С, является гармонической и 
теряет это свойство в пространстве х, у, і, если 
Я,^Я 2 ^Я 3 (Ѵ 2 р{і;,гі,ф=0, Ѵ 2 р(х,у,і)ф 0). Из этого сле- 
дует, что система (5) при ѵ(%,т 7,4), р(%,гі,ф не отра- 
жает законы гидромеханики в обычном понима- 
нии, и ее возможные решения при данной зависи- 
мости дают лишь приблизительную информацию о 
реальном явлении. 


Обтекание эллипсоида. Скоростное поле 
и поле давления 

Если эллипсоид находится под^ углом атаки, то 
скорость набегающего потока Ѵ(\Ѵ\=\) удобно 
представить в виде разложения: 

ѵ = ѵ 1 і + ѵ 2 ]+ ѵ і іс , е; 2 +к 2 2 + е 3 2 = і. 

С тем, чтобы удовлетворить уравнению несжи- 
маемости, скоростное поле ищется в виде: 

Е = = ~гоТ(/Ехсг), / = /(а), сг 2 =| 2 +ц 2 +С 2 , 

а( 2 

где /(&) выбирается из требования выполнения 
последнего уравнения системы (5). 

Расшифровка операторов пространственного 
дифференцирования показывает, что 

ѵ=(/ + ^-а/^Е-|(<7К)а~ 1 /'<7, 


ГОІ ГОІ ѵ = - 


гоіѵ =— — (/' + 4<т '/')Ѵу.6, 


(а(Г + 4а-'Г)' + 2(Г + 4а-'Г))Ѵ- 
-(<7 Ѵ)а' (/” + 4а'/')’а 


ГОІ ГОІ го! V = —[(/" + 4(7 1 /У +4(7 1 (/"+4сГ ! /У]КХ(Т, 
гоі(ѵхгоі ѵ) = -^(дѴ)а-' /(/" +4а-' Г)'ѵ ха. 


Теперь видно, что третье условие (5), соответ- 
ствующее требованию (3), будет выполнено, если 
функция /ст) удовлетворяет обыкновенному диф- 
ференциальному уравнению: 

(/" + 4а- 1 /')" + 4а- 1 (Г + 4а- = 0 , 


что в общем случае приводит ее к представлению 
по (4). Заодно легко просматриваются варианты 
для исполнения полного уравнения Гельмгольца, 
получаемого после применения к первому соотно- 
шению из (1) операции ротирования. 

Если при внешнем обтекании потребовать, что- 
бы при о=1 поверхность объекта (в пространстве х, 
у, і им служит эллипсоид, а в пространстве г), / - 
сфера) была поверхностью тока и на ней осущес- 
твлялось прилипание, а вдали от него было 

-/ ч т? Г- 4а 4г 

ѵ(со) = V ѵ = — ф — 

V ей й( 


то (4) примет вид: 

/' = 1-— а -1 
2 


1 


—а 

2 


а это приводит к тому, что скоростное поле предс- 
тавится равенствами: 

-^К 2 |7 7 а- 3 (1-а- 2 )-^Ез^7 7 а- , (1-а- 2 ), 

Ц =-^<5ца- 3 (1-а- 2 ) + 

+К 2 (і-^а-'-іа- 3 -^ 2 а- 3 (1-а- 2 )^- (6) 

^=-^^а- 3 (1 -а- 2 )-^ 3?7 Са- 3 (1-а- 2 ) + 

+ Г( 1 - 

А поскольку теперь 

гоігоі ѵ = ^(а" 3 Е-3(а Т)а" 5 а) = ^Ѵ((а Е) а -3 ), 

то из закона движения по Стоксу поле давления бу- 
дет определено формулой 

р = -^(аѴ)а\ (7) 
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По форме система равенств (6) и (7) выглядит 
совершенно одинаково как для эллипсоида, так и 
для сферы. Но в первом случае надо читать: 



а во втором: 



У_ 
К ’ 


с = 


к’ 


где а, Ъ, с, К- соответственно полуоси эллипсоида 
и радиус сферы. Вполне понятно, что при наличии 
только одной составляющей скорости набегающе- 
го потока, направленной вдоль оси эллипсоида, 
формулы упрощаются, а для сферы они воспроиз- 
водят известное [ 1 ] решение Стокса. 

Еще раз следует отметить, что найденные поля 
(6) и (7) для обтекания эллипсоида не являются 
вполне достоверными, так как безразмерная систе- 
ма (5) при учете масштабов а, Ь, с будет выглядеть 
несколько по иному. Вопрос о степени достовер- 
ности найденного решения можно выяснить после 
вычисления сопротивления взятых объектов. 


Сопротивление сферы и эллипсоида 

Сопротивление состоит, как известно, из двух 
слагаемых - сопротивления от давления Р(р)=-\рсіР 
и сопротивления трения Ну)=2\усІ$ . Здесь через 
обозначен направленный по нормали поверхност- 
ный элемент объекта, а через у - тензор скорости 
деформации, компоненты которого в силу отме- 
ченных особенностей скоростного поля выглядят 
совершенно одинаково для сферы и эллипсоида. 
На самой поверхности (при ст=1) они запишутся в 
виде формул: 

г«= 

У^=\(Ѵ,ч(\-2?) + ѴЛ(\-2ц г )-2Ѵ ъ $ц(;Х 

/ к =|(^С(1-2Г)-2С^77С+С 3 |(1-2С 2 )), 

Ущ =|'7(-Е;<5і7 + С 2 (1-ц 2 )-Сз77СХ 

У га = | (~2К Ж + Ѵ 2 С (1 ■ - 2г ] 2 ) + Ѵ 3 77 (1 - 2С 2 )), 

г«=%С(-г,%-ѵ 2 т 1 с+ѵ 3 ( 1-С 2 ))- 

При этом давление на поверхности принимает 
значение: 

Р=- |<тК = -|(С^ + С 2 т 7 + С 3 0. 

При производстве процедур по интегрирова- 
нию необходимо учесть, что 

| = СО 5 0 , Г] = 8 ІП Ѳ С 08 ф, С, = 8 ІП 08 ІП(р; 0 < Ѳ < Л, 0 < (р< 2 п, 

благодаря чему поверхностные элементы берутся в 
виде записей для сферы и эллипсоида соответственно 

сІБ = (<* / + Ц У + <= к ) 8Іп Ѳ сІѲ сІ(р, 
сіЗ =(Я 2 Я 3 §/ +Я 1 Я 3 ?7у +Я 1 Я 2 4’Я)8 ІпѲ сіѲ сіц>. 


Исполнение процедур по вычислению приво- 
дит к результатам для сферы 
Р(р) = 2лѴ, Р(у) = 4лѴ, Р{р,у) = Р(р)+Р(у) =6 лѴ, (ІЙ =1) 

и для находящегося под углом атаки эллипсоида 
Р{р) = 27г(Я 2 Я 3 Ѵ х і + Я 1 Я 3 Ѵ 2 ] + \Х 2 Ѵ 3 к), 

Р(у) = |я:[(4Я 2 Я 3 +ЗЯ,(Я 2 +Я 3 ))И/ +...] (по циклу), 

Р(р,у) = -|л-[(ЗЯ 2 Я 3 +Я,(Я 2 + А 3 ))Й|/ +...] (по циклу). 

Итоговая информация в размерных величинах 
получается после умножения соответств^?ющих_вы- 
ражений на размерные комплексы ц| Ѵ.,\Р, /и\К\Ь, 
благодаря чему она фиксируется формулами: 

Р(р,у) = 6лрКѴ хі , 

(ЗЪс + а(Ъ + с))Ѵ 1х і + 

+(3 ас + Ъ(а+с))Ѵ 2х ] + 

+(3 аЪ + с(а + Ъ))Ѵ 3а> к 

соответственно для сферы и эллипсоида. Для вари- 
анта обтекания без угла атаки, например, при 
Ѵ. = Ѵ„і итоговые формулы понятным образом уп- 
рощаются и принимают вид: 

X* = 6лрКѴ„,Х, = |тг ^Я(3 Ъс + а(Ь+с)). (8) 

Обычно сравнение сведений по сопротивлению 
производится путем введения эквивалентной сферы, 
радиус которой К=К(а,Ъ,с,Ь ) получается в результате 
совпадения сил X, и X' из формул (8), что дает 

1К а =РрЪс + а(Ь + с)1 

Индекс подчеркивает, в направлении какой из 
главных центральных осей инерции движется эл- 
липсоид. Здесь масштаб обезразмеривания I остает- 
ся неопределенным, и способы его выбора нуждают- 
ся в дополнительных обоснованиях. Наиболее прос- 
тым (напрашивающимся) выбором этого линейного 
размера является равенство Ь=К, что влечет за собой 

К 2 а = 1[3 Ъс+а(Ъ+с)1 ( 9 ) 


И 


здесь и далее предполагается, что полуось парал- 
лельна направлению набегающего потока. 

Но следует заметить, что такой выбор размера Ь 
в случаях вырождения эллипсоида в пластины, ког- 
да преобразование координат перестает работать, 
вряд ли можно считать оправданным. Например, по 
теории потенциала притяжения известно [1], что 
предельные выражения для сил сопротивления 
плоского круглого диска (радиуса й), движущегося в 
направлениях, перпендикулярном своей плоскости 
и параллельном ей, соответственно равны 

Х = 16 рЪѴ, Х = ^-/лЪѴ, 


что при сравнении с эквивалентной сферой означает 


К 

~Ь 


— я 0,85, - 
Зтг Ъ 


фб 

9л 


0,57. 
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Рисунок. Отношение эффективного радиуса сопротивления к величине полуоси: а) поперечной направлению движения, б) па- 
раллельной направлению движения 


В то время как из (9), в случае й=0 и Ъ=с, соот- 
ветствующем варианту движения диска перпенди- 
кулярно своей плоскости, имеем 


К 

ъ 


л/15 

5 


0,7746, 


а при а=Ъ и с=0, что соответствует движению па- 
раллельно плоскости диска: 

п Гс 

— = = 0,4472. 

Ъ 5 


Попутно можно отметить, что по [1] сопротив- 
ление диска при движении в ортогональном к сво- 
ей плоскости направлении в полтора раза больше, 
чем в продольном, а по (9) это отношение состав- 
ляет Ѵ3=1, 732. 

Подобно [3] для течения, параллельного оси 
симметрии (при Ъ=с), введем отношение длины к 
диаметру е=а/с. Из (9) отношение эффективного 
радиуса к величине поперечной полуоси будет 

А = ^=^(3 + 2е). ( 10 ) 


Из [3] данная величина при е>\ равна: 


С- 6 

3 


1 


2 Е 

Б 1 - 1 


2е -1 


(е 2 -1) 


, 3/2 


1п 


Б + 


4б^і 




а при е<1: 

М 


1 


2 Б 

1-е 2 


2(1 -2е 2 ) 
(1-е 2 ) 3/2 


агс!§ 


л/^ 


(П) 


На рисунке, а приведено сравнение формулы (9) 
с данными [3] при различных значениях е при дви- 
жении эллипсоида вдоль оси симметрии. В случае 
движения поперек направления оси симметрии 
следует положить а=Ъ и ввести е=с/а. Тогда (9) дает 


К = 


К 


(4 е + 1). 


( 12 ) 


Согласно [3] при е>1: 

О 

5 = 

2е 2 - 


е -1 (е -1) 


угу 1п(е + л/е 2 - 1) 


а при е<1: 


5 = 


8 

3 


1 




(13) 


На рисунке, б приведено сравнение при различ- 
ных значениях е в случае движения тела поперек 
оси симметрии. Как видно, (9) дает всюду зани- 
женные результаты при хорошем совпадении для 
тел, мало отличающихся от сферических. Следует 
отметить, что при е— >0 формулы (10) и (11) соотве- 
тствуют диску, движущемуся в потоке, перпендику- 
лярном своей плоскости, а формулы (12) и (13) - 
диску, движущемуся параллельно своей плоскости. 

Полученная информация о сопротивлении при 
значительных отличиях полуосей а, Ъ, с от радиуса 
сферы К, существенно расходится с имеющимися 
теоретическими сведениями. Объясняется это, по- 
видимому, тем, что использованная система (5), за- 
писанная с применением изотропного обезразме- 
ривания, некоторым образом отличается от той, 
которая получается при неизотропном обезразме- 
ривании. Это обстоятельство приводит к тому, что 
введенная функция Дет) должна иным образом за- 
висеть от аргументов г /, С,. 


Заключение 

Таким образом, решена задача о медленном об- 
текании эллипсоида, мало отличающегося от сфе- 
ры, указана формула для расчета его сопротивления. 
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Использованные обозначения 

ѵ -вектор скоростей, Ей - число Эйлера, Ее - чис- 
ло Рейнольдса, р - давление, у/ - функция тока, 
Ѳ - сферическая координата, ср - потенциал, а, Ъ, 
с - полуоси эллипса, т). С, - криволинейные ко- 
ординаты, х, у, і - декартовы координаты, /и - ди- 
намическая вязкость, X, X, - силы сопротивления, 
с!5 - элемент поверхности. 
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